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The perturbative approach for describing the underdense plasma-ultrarelativistic elec- 
tron bunch system is developed, using the ratio ^ as a small parameter (n5-bunch,r2o- 
plasma electron densities). Focusing of the electron bunch emerged in the first approx- 
imation of the perturbative procedure as a result of the plasma electrons redistribu- 
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' tion. Focusing gradient and strength for ultrarelativistic,fiat, uniform and short bunch are 

O '■ 

\0 • obtained and compared with the previous results. 
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i INTRODUCTION 

Plasma focusing devices, being compact, simple and effective elements are promissing for 
obtaining electron (positron) beams of very small spot size, requested for future high 
■ energy linear colliders. Theoretical predictions and investigations of the plasma lenses in 

overdense {ub < uq = Up) and underdense {ub > nQ,nb-hunch electron density,no-cold 
neitral plasma electron density), performed last decade |T[|-[|| and were supplemented by 
experimental tests, carried out at ANL 0]-|Ty],Tokio University-KEK [ll|] and UCLA [12 



for overdense plasma lens regime. 

The theoretical treatment of the overdense plasma lens,performed in linear approxi- 
mation,using ratio ^ as a small parameter, is more or less complete, at least in the frame 
of the cold plasma and rigid electron bunch approximations. 

The existing theoretical approaches to underdense plasma lens are phenomenological 
by nature. 
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It seems, that more detailed description of the underdense plasma lens is needed, in 
particular, taking into account the program for the future experimental investigations 
[l^, which include underdense regime too. 

Using the ratio — <^ 1 as a small parameter, and developing the subsequent per- 
turbative approach for the description of the underdense plasma lens, it is possible to 
achieve, at least the same level of the understanding of the focusing phenomenon in un- 
derdense regime as that of overdense plasma lens case. 

The present work is devoted to this goal. In the next section the perturbative approach 
is developed for underdense plasma regime for the system of the Maxwell equations and 
hydrodinamical equations of the motion of the plasma electrons (ions are assumed im- 
mobile). The third section devoted to zero approximation calculations of the generated 
fields and plasma electron velocities, assuming that bunch is ultrarelativistic and short 
enough. In the section 4 the first order approximation results for plasma electrons redistri- 
bution and transverse force (sum of the electric and magnetic forces) are given. Section 5 
contents the dissension of the obtained results. In particular, it is mentioned that physical 
picture of focusing phenomenon in underdense case is based on the redistribution of the 
plasma electrons around the driving bunch. The effect of the ion column in ultrarelativistic 
case is negligible, compared to the focusing due to redistributed plasma electrons. 

The mathematical technique developed and above mentioned physical discription of 
the focusing phenomenon may be used for more detailed, coherent and associated with 



the experimental program |13] computer calculations 



A PERTURBATIVE APPROACH FOR UNDERDENSE PLASMA LENS 

For definiteness, consider the flat electron bunch with the vertical dimension 26,which is 
assumed much smaller than horizontal dimension 2a; longitudinal dimensions,which are 
arbitrary,are 2(i.The bunches of such a geometry are suitable for future high energy linear 
electron-positron colliders |jl3[ . 



Bunch electrons uniform density is and bunch is considered as ultrarelativistic 
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and rigid one. Plasma electrons density is uq <^ nfe,plasma is neutral,cold and with the 
immobile ions. 

The geometry of the bunch,which moves in lab system through the plasma with 
constant velocity fo, allows to consider the electric and magnetic fields components as 
follows: 

= 0, Ey, E, ^0;B,^0,By = B, = 
This more or less evident approximation for the fiat beams with large aspect ra- 
tio, along with the condition -i-. ^ -i- was introduced in [0.1n |TH] it was pointed out, 



dx dy 

that the assumption By <^ = B "is not universally true everywhere inside the even 
very fiat beam" .It is neccessary to take this into account at precise analytical or computer 
calculations. 

All physical quantities of the problem are considered as a functions of the argu- 
ments y and z = z — v^t only (steady state regime). Introduce the dimensionless ar- 
guments y' iz' = kby,kbz;t' = ooht^uol = , kb = ^,and dimensionless variables 
E'y,, = i^y'Eyy.B' = i^r'B.,; 

Followig ||I6|, |17|, define the generalized plasma electron velocities as 



T/ _ ^ey _ 13 ey ft = US. fi — ^e.zi; T A _ V^z _ Pez 

vo-vez P-Pe.'^ c^^'^^^y C ' VO-Vez I3-I3ez 

and generalized plasma electron density as 

(BTFCh transformations |]ro|),it is possible to rewrite the continuity equation in the 
following form: 

-f + ™^0 (1) 

oz oy 

(in (|1]) and what follows the tilda and prime superscripts are removed). The Maxwell 
equations in this case are: 
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1. 

2. 
3. 
4. 



dB 

dy 
d 



(5 + (3NV, + /5 



(92; ' 



^iB + PEy) = -l3NVy 



(2) 



+ 



[l-noM + iV(l + V;)], 



^2; 9?/ 

Eqs. (]^) are valid for — d < z < d, —b < y < 6;when \z\ > d, \y\ > b in right hand side 
of eqs. (0.1) and (@.4) the quantities /? and —1 subsequently are absent. Eq. (||.4) follows 
from ®, (0.1),(0.2). 

The eqs. of the motion are: 





Vz = 




1. 






Vy = 




2. 




where 







-E,{1 + 2\4 + ^) + PXV, (^Ey + ^B]+ (3VyB 



(3) 



V=1 + 2V, 



2t^2 



(4) 



^ dz ^ dy ' ^ dz 



y ' 



dy 



(5) 



Adopting the condition — <^ 1, decompose the quantities in question in the following 



series: 



N = eNi + e^N2 + ... (6) 
Vz = Vzo + eV,i + e^Vz2 + ... 
Vy = Vyo + eVyi + e%2 + ... 
Ey^z = EyzQ + ^Eyzi + e^Ey^z2 + ■■■ 
B = BQ+eBi+e^B2+ ... 
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where e = ^ (in what follows e" is included in subsequent quantities A^„, Vm and so 

on) 

In the zero order approximation from (H.2) and (0.3) 



Bo = -PEyo, = (7) 

From (|^ and (0.4), introducing the potential ipo{z,y) by 



- Eyo - (8) 



it follows that 



1 aVo 5Vo \ h -d<z<d, -l<y<b 

+ ^ = (9) 



7' dz^ dy^ [ 0, \z\ > d, \y\ < b 
In the first order approximation: 



dz dy 
from (lOD and (|) follows: 



+ l^r^ = T-[ 7^1-^iKo ; (11) 



The focusing force is 



/, = -e{Ey + (3B) = fyo + /,! + ..., (12) 
fyo = -e{Eyo + (3Bq) = -e'j'^Eyo, 
fy, = -eiEy,+pB,) = -eWyi. 



From (0.3) in the first approximation 



dWyi ^ dE^ 
dz dy 



(13) 
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which is an analog of the Panoffsky-Wenzel relation;remember that in ( p!3[ ) z = z — v^t. 
Different ating (0) over y and then integrating it over taking into account relation (|l^) 
and that 



A ~ ^^^1 ~ A^i^^o) when y ±oo,we have 

^J^A.L^^(^-,^N.-N.V^ (14) 

where z = ^z. 

In what follow,we consider the ultrarelativistic bunches,when 7^1; then /^o ~ 



0(7 ^) and in right hand side of ([I^) it is possible to leave only the term — So for 



(9y 



the sought quantity Wy\ we need A?"! from eq. (p^Ol). Entered in (|TOD Vy^ must be found 
from eqs. (^,(^ for the fields in zero approximation.lt is evident,from eq. (|T^,that the 
role of the noncompensated positiv ions (the term on the right hand side of eq. (0), 



proportional to ^ — A''i(l + Vzo)),is negligible in ultrarelativistic limit (~ 0(7 )). 
3 ZERO ORDER APPROXIMATION 

In the zero order approximation it is necessary to solve eq. (|^) for the flat bunch moving 
in vacuum {J^ = 0). 



The solution may be written using Green function formalism (see e.g. |T8 



1 r<^ 



^,{z,y) = ^ I _dz' I dy'lnUz-zf + {y-y'fr^ (15) 
Zir J-d J-b L J 

where d = •yd. 

Double integration in ([l5| ) may be performed using standart integral tables, but the 
result of integration has a complex form,which is difficult to use. It is more convinient 
to use (^ and perform an approximate integration for 7 ^ 1, taking into account, that 
argument of the integral has a complex pole, when —d < z < d, —b < y < 6. Under these 
conditions the integration over z along the real axis in complex z plane is equal to residue 



1/2 
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at the pole, minus integral over the semicircle with the radii d in upper complex z half 
plane. The approximate results of integration for 7 ^ 1 are the following: 



—b<y<b, —d < z < d; 



y >b, —d < z < d; 



E,o - (18) 
' Try 

y > jd b, —d < z < d; 



E,. « (19) 



z:^ d,'jd:$> b,-b <y <b] 



E.0 = - In \y-bf « l'\z - d\' (20) 

7r7 (a + z) 

E,o = -^2^, \y - hf ^f\z- y>-id^b 

T^y 

For comparison, remember that for the point relativistic charge 



Ey = E^ = ^7, E, = E\\ = ^,R^ = x' + y' + {z- Votf; (21) 

i.e. the electric field for 7 ^ 1 is practically transversal. In our case of the flat 
bunch expression (p^pOD shows that always \Ezo\ <^ ji^yol^but the expressions for Eyo, E^q 
modified, compared to ([21|),due to the different geometry of the charge distribution. In 
what follows, we will use the approximate expressions for Eyo up to terms ~ 0(7"^), and 
put it to zero when \z\ > d. 

The next problem in zero approximation is the definition of the generalized velocities 
Vyo, Vzo-When 7 ^ l,and bunch length is short enough d' <^ 1, {d' = kbd) it is possible 
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to drop out in the terms proportional to \Vzof <^ \Vzo\ and to \Vyof ^ |K/o|- Taking 
into account that in the considered case \Eyo\ ^ \Ezo\ the equations of the motion (||) 
approximately can be rewritten in the following simple form 



^ - VyoEy,, (22) 

dV 

-g^^il + V,o)Eyo,\y\<b,\z\<d; 
The solutions of eqs. (^) for the boundary condition Vyo = 0, Vzo = 0, when z = d 

are: 



Vyo = sinh {\Eyo\{d - z)) ~ \Eyo\{d~z), (23) 
Vzo = cosh(|^j^o|(c^- -2)) - 1 ^ -\Eyof{d-zf, 



where EyQ is given by rtT6-18)and Vyo — > 0, V^o — ^ when y ±oo. For the short 



{d ^ 1) bunches from (^) follows that V^o ^ VyQ. 
4 FIRST ORDER APPROXIMATION 



Using eq. (|l]) in the first approximation (see ([T0|) ) and eq. (p3|) for l^o it is possible to 
find Ni{y, z) by the method of characteristics. 

When \y\ < b, EyQ ~ — y + 0(7^^) and characteristics are 

ye-W~-f = C,,yNi'^ = C2. (24) 

The solution of eq. (p^), which is equal to when z = d,is 
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When |?/| > b, EyQ ~ — 6 + 0(7"^) and characteristic is —^{d — z)"^ + f = C,so n[^^ is 
function oi —^{d — zf + j^; When y = b it must coincides with the expression (P5D,so 
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(2) = I^e-^^-^-^^'+^-i 



rib 



When y > "yd ^ b, Eyo ^ — and characteristics are y'^ — ^^{d 
Taking into account that N^^^ — > ^ when \y\ ^ oo,we have 



(3) 



(26) 



Co 



(3) _ '^O 



\y\ 



Tib 



1/2 



(27) 



y2-?>^^d~z) 

Expressions ( P5[p7] )are the solutions of the eq. ( pij ) for subsequent values of Vyo 
^). Expression (|26D for z = d gives 



(28) 



which seems as unappropriate. Taking into account, that (|26D is an approximate expres- 
sion,we assume the following interpolation formula for . 



(2) 



Ub 



(29) 



which is — when z = (i,and coincides with the 

rib ' 

(2) 

Ni is exponentially rising,when y — 



when y = b. 



+00 and at the some point ycr it intercept 
the Ai'a^'' , which is also rising,when y yp + = ^ {d — z) < ycr- Hence A^i has a 

maximum value and the plasma electrons distribution exibits a crest at the point y = ycr- 
Parameters c and ycr may be defined by the condition 



Ni'\ycr,z) = Nl'>{ycr,z) 



(3), 



(30) 



The second condition may be obtained,using the following consideration: longitudinal 
displacements of the plasma electrons are negligible, compared to transverse displace- 
ments, due to condition Vzo <^ Vyo,so for the given \z\ < d deficiency of the plasma 
electrons in the region < y < ycr (and in the symmetric region — < y < 0) is equal 
to the surplus of the plasma electrons in the region ycr < y < +oo (and in the region 
—oo<y< —ycr subsequently): 
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!^-0<i,+ rf^-A'P'H!/= (31) 



\nb J Jb \nb 

Jyc J Jycr V fib J 



In (^Tj) yc is the point, where iV(2,^) = ^ From for d < 1: 



(1 



„2 



From eq. (pTj) for d <^1: 



c=- \ (32) 



= ^ » 6 (33) 

TT 



and c from (^) approximately has a following simple form: 



Now it is possible to turn to determination of the focusing force Wyi from eq. (]T^). 
For 7 ^ 1 eq. (0) we get an approximate form 

dz'^ dy"^ dy 
and using the Green function formalism the solution of eq. (|35|)can be written as 

Wyi = -— dz dy — \n [y -y) + [z - z) . (36) 



27ry 7 ^ ay' 

The domain of the integration in (^) is defined by the condition ^^gj'^"* 7^ 0,i.e. 
it is —00 < y < +CXD, — 00 < z < (i,because for z > d, Ni = ^ and = O.For 
—d < z < d,^^ is given by p5|),([2^),(p^);for z < —d from the continuity condition 
Ni{z, y) = Ni{—d, y) up to some —zq < —d where considered steady state regime changes 
to (nonlinear) wake wave (and later on the uniform distribution of the plasma electrons 
with A^i = ^). 

^ rib' 
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Due to the symmetry of the Green-function and antisymmetry of the derivative 
it is possible to change the integration over the region — oo < ?/ < to the integration 
over the region < y < oo. Then the expression (|36| ) has the following form: 



Alt 1 -In h " dy' (y + ylf + (=-z'f ^ ' 

taking into account (|25|).The interval of the integration over y' in ( p7D must be divided in 
the two intervals (6, ?/cr) and (?/cr, oo\ where eq. (^) and (|27|) must be used for ^^;both 
expressions (p7| , [29|) have a maximum at y' = ycr and it is possible approximately to 
take out of the integral the slowly variing logarithmic function at the maximum point 
y' = l/cr- Then (|37|) will have the following form 



Tib 



(38) 



Notice that the approximate expression ( p8| ) for Wyi does not depend exphcitely on the 
form of, then to some extent arbitrary,function N^^^ (pOf), because N^^\h,z) = A''}^^(z),see 
(^) and (1^). Only the value of ycr (|33[) , which enters in (p8|), explicitely depends on the 
form of Ni through conditions ( pO| , pl|). 

Taking into account, that \y\ < b ycr it is possible to develop the logarithm function 
under integral in a Taylor series, leaving the first term of the expression; then for ci -C 1, ?/ < 
b the expression for Wyi will have the following form: 

1 ^0 , > id- z'f , , 

27r nfo 7-^0 yl^ + (-2 - z') 7^ 
Consider first the integral over z' in the domain (— d, (i),and later on in the domain 
(— ^0, — (i).The resulting expressions denote by and ^ subsequently;then 

(0) ^ 2a(^!^^^2^ (40) 

where q;(z) is known function, which varies in interval 0,21 < ai^z) < 2,26,when 
—d < z < d,Wyi^ ~ ^(3{z, zo)^yd'^, (3{z, Zq) weekly depends on the value of the zo^and 
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it is possible to take for estimates zq +00. Then (3{z) = (3{z,oo) varies in interval 
0, 26 < (3{z) < 1, 57,when -d < z < d. 
Hence, for estimates, it is possible to use 

Wy, = W^? + Wi^^ ^ -yd'[a{z) + P{z)]^ ^ -yd^z)- (41) 
" ^ -n rih IT rih 

It is possible to apply the same approach to the flat, rigid, short positron bunch with 
uniform charge distribution. In this case in eqs. (|15[|T^) electric field components change 
sign. As a consequence the Voy component of the generalized velocity (eq. |2^) also changes 
sign,which corresponds to attraction of the plasma electrons by positron bunch. Obviously 
enough, this leads to drastic change in plasma electrons redistribution as compared to 
electron bunch case. Instead of eq. (|25|) for \y\ < 6 it is possible to obtain 



Nff = !^e^('^-^)^ (42) 
rib 

instead of (|2^) ioi y >h we have 

ATg) = "^^W-^f+y/^-^ (43) 
rib 

and instead of ( pT] ) for y > ■yd b 

< = ^ ^ n,; (44) 



From (^) , (^3|) it is seen that plasma electrons concentrate around the rear part of 



the positron bunch. From ( ^4]) it is seen that N^p' has a complex pole and hence is slowly 
rising function, approaches it's limit equal from bellow,when y — > +00. This means 

Til) 

that plasma electrons came to the region occupied by positron bunch, \y\ < 6,from the 

regions,where y > yd ^ 6. Focusing effect for the positron bunch case is caused also 

(2) (3) 

by plasma electrons redistribution. But due to different behaviour of N^p and N^p , the 
interpolation between these two,which provide the main contribution to focusing force 
(p!^),(^7D,is much more arbitrary,than in the electron bunch case. In any choose of the 
interpolating function it must be more or less smooth, and therefore the estimates of the 
integrals like (|36| ) , (|3^) will be cumbersome enough. 
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These two obstacles forced to think that positron bunch case must be treated numer- 
ically,even in the short bunch case and in the presented formahsm. 



5 DISCUSSION 

The expressions for the focusing force, focusing gradient and stregth are subsequently 

fly = -eW^y, G = ^,K= (45) 

Wiy is given by (^T]), which is valid when 7^1 and d' -C l{kbd <^ l).In the usual 
units expression (|4l|)has the form 

Wy, = ^^^^{non,){yd') (46) 

TT 

where = e'^/mc^ is the electron classical radius. Subsequently,focusing gradient and 
strength are: 

G = ^^(n,r^o)c^^ K = ^^^(n^no) (47) 

TT 77r 

The expressions {^7\ } differs from the previous results,which for the flat, long enough 
electron bunch and underdense lens are 

G = A7rnoe,K=^^^ (48) 

7 

These are based on the understanding the focusing phenomenon as caused by the positive 
charge of the ion column, which exists inside and behind the moving electron bunch, 
because the considered bunch blows out the plasma electrons. 



As it follows from the presented consideration (see,e.g.,eqs. (14, ^) the net focusing 
effect of the ion column, caused by noncompensated positiv charge ^ — Ni(l + Vzo), is 
proportional to 7"^ and is negligible in our approximation 7 ^ 1, compared to focusing 
effect of the plasma electrons redistribution, described by eqs. (^) and (|39D .The presence 
of the factor 7"^ in the part of the net (electric plus magnetic) focusing force, which caused 
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by ions, is a result of magnetic compesation of the subsequent focusing electric field in lab 
frame, where electron bunch moves with the velocity Vq and unperturbed parts of plasma 
are at rest. Indeed, from Maxwell eqs. (0.1), (0-4) in the first approximation: 



dB 



dE. 



■zl 



(49) 



dE. 



■zl 



+ 



dy 
dEy, 



iVi(l + Ko), 



dz 



dz 



dy 



it follows: 




entered in (^91), represents the nonuniform noncompensated positiv ions current, which 
propogate with the velocity fo,equal to the velocity of the electron bunch. This term 
caused the compensating magnetic field,which at 7 ^ 1 cancels the focusing action of 
the subsequent electric field. The physical reason for that - the exposed positiv ions col- 
umn effectively "moves" in plasma along with the generating electron bunch with the 
velocity of it. When bunch, moving through plasma,blows out the plasma electrons, it un- 
cower new forward parts of the partialy noncompensated ion column and this "motion" of 
the revealed positiv ions charge, i.e. subsequent current, generated compensating magnetic 
field. 

The physical interpretation can be based also directly on the displacement cur- 
rent, entered in the first eq. (^) through the term /3^^.Ions ahead of the bunch are 
completely compensated by plasma electrons, inside the bunch ions charge is only partially 
compensated (eq. (p5l)),so exists the change of electric field component with z = z — vqI, 
which generates subsequent compensating magnetic field. Formaly this part of displace- 
ment current coincides with the current conditioned by "effective ion motion", as it can 
be seen by replacing the electric field derivative in the first eq. (^) by it's value from 
Coulomb low-second eq. (|49|) .It is evident that this coincidence is based on the steady 
state regime, adopted at the present work. 
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In the bunch frame magnetic forces acting on the bunch electrons are absent. Electric 
field of the bunch charge has no relativistic contraction and hence the plasma electron 
redistribution has a different shape,more uniform,than in lab frame,which dimenish their 
focusing effect. But the focusing effect of the plasma electrons still is proportional to 
rife-bunch charge density;their focusing electric field is relativistically contracted. These 
factors increase the redistributed plasma electrons focusing effect, which will be larger 
than focusing effect of the electric field of the ion column, proportional to no <^ nfe.Lorentz 
transformation did not changes the transverse component of the focusing force, so it must 
be the same in both frames. Above mentioned physical arguments allow to think that it 
could be the case. It is difficult to perform analytical calculations in the bunch frame due 
to complexity of relativistic equation of motion of the plasma electrons. Recall, that in 
presented calculations the plasma electrons considered as nonrelativistic,due to condition 

(M) < 1. 

All above mentioned general physical arguments are valid for the relativistic bunches 
with arbitrary length and charge distribution. Assumptions adopted in the present work 
(short enough relativistic bunch with the uniform charge distribution) allows us to per- 
form analytical calculations. 

Unfortunately the condition kbd ^ 1 (rf' ^ 1, short bunches) does not permit to use 
expression (^Tf ) for focusing strength for the flat beams of the FFTB [|13|.In that case 
rife = 7, 7 ■ 10^^cm~3 and 2, 8 ■ lO^^cm-^ and h = 5,2 ■ 10^cm~^ and h = 3,1 ■ lO^cm-^ 
consequently, d = 2,3- lO^^cm and hd = 1, 2 ■ 10^ > 1 and 0, 7 ■ 10'' > 1. 

Using the short bunches, nevertheless it is possible to obtain sufficiently strong focusing 
gradients. For example let k^d will be kbd = 0, l.Then d = 10^ ( ) ,G = 7, 5-10^ no 



and if nn ~ 10^'' cm~^ as in FFTB case,G ~ 7, 5^^. In order to fullfil the condition 

" ' ' cm 

^ <^ l,nb must be e.g., equal 10^^cm~^,then d = 5 ■ 10~^cm.May be it is possible to 
obtain such a short bunches experiment ally, using the technique, proposed in [^,|^. 

It is worthwhile to stress, that due to specific form of the redistribution of the plasma 
electrons focusing gradient ~ and focusing strenght ~ have a quadrupole charac- 
ter, due to nonuniform shape of the distribution of plasma electrons along the crest, which 
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diminish, inspite of the factor (nonf,), significantly the numerical values of these quantities. 

This disadvantage probably can be avoided, considering the medium kpd > 1, or long 
kpd ^ 1, bunches. The main difficulty in the analytic approach to these cases lies in 
the complex form of the equations of the motion |^) which are fully relativistic in 
these cases. Probably the medium and long bunch cases can be treated by computational 
methods only. Presented analytical approximate description of the short bunch case can 
be used then as a physical guide and programs testing example. 

It will be interesting to test experimentaly the predictions ( ^7D for focusing gradient 
and strenght,particularly,the dependense of these quantities on rih and d. 
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